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Abstract 

In order that current and future renormalon results in QCD can 
be used to their fuh advantage it is important to understand how the 
Borel transforms of related functions are themselves related. For ex- 
ample, a change of renormalisation scheme can involve the inversion 
of a function; how is the transform of a function related to the trans- 
form of an inverse ? We outline a method of extracting the dominant 
behaviour of the Borel transform after the action of simple operations 
on the function. 
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1 Introduction 



Borel transforms provide a convenient and insightful method of describing 
the large-order behaviour of QCD perturbation series. An observable 



ra=0 



R{a) = J2 rnaJ'^\ a = a./vr (1) 

is represented as the integral 

R{a) = I F{z)e-'/''dz (2) 



where F{z) is the transform. (If necessary, principal values are used to define 
the integral - see Appendix). A singularity of the form 

gives rise to a factorial divergence 

nl 

in the perturbation expansion of R{a). 

Much attention has been devoted to characterising the F{z) of QCD 
observables, both in general and for particular cases. Arguments related 
to power corrections show that a typical observable has an infinite chain of 
singularities of the form equally spaced along the positive real axis of 
the Borel z-plane [00. These are called IR renormalons and are crucial 
is assessing the practical reliability of perturbation theory. Although the 
positions and general form of the renormalons are fixed for any R{a) by these 
arguments, the particular details (especially the A and /5 of a singularity) 
have to be established for each observable separately. Such work has both 
reproduced the predicted general features and provided specific information 
about F{z) in cases of interest, see e.g. [^-0. A notable feature of these 
results is that F{z) is invariably a simpler function than R{a). 

However to fully exploit the potential of these results it is often necessary 
to consider functions related to R{a) rather than the calculated R{a) itself. 
For example, R{a) and its Borel transform may have been calculated using 
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the MS renormalisation scheme, but it is now the corresponding resuhs in 
the effective charge scheme [10| that are wanted. The function of interest is 
now the effective charge /3-function 



p{R) 




(5) 



where f3{a) is the MS /^-function. Calculating this involves doing a series 
of operations to the function R{a): inverting it to give a(i?), finding the 
function-of-a-function f]{a{R)), differentiation and division. In terms of the 
function these are straightforward enough operations, though possibly diffi- 
cult to do exactly in practice. A natural question - particularly if, as noted, 
the observable is simpler transform than function of a - is what effect 
these operations have in the Borel plane ? If we know the Borel transform 
of R{a), what does this tell us about the transform of a{R) ? 

Because such operations are commonplace, questions hke this arise con- 
tinually when applying the results of renormalon calculations. It is this 
practical need that the method outlined here seeks to meet. 

2 The Operations 

The operations that this paper will discuss are 

1. Multiplication Ri{a)R2{a) 

2. Differentiation dR/da 

3. Function-of- Function i?i(i?2(o)) 

4. Division i?i (a) /i?2 (a) 

5. Inversion a{R) 

Some of these are easy and some are hard. The first two are the straightfor- 
ward ones. Like with any simple integral transform, the Borel transform of a 
product is given by a convolution integral over the transforms of the factors. 
If Ri{a) and R2{a) have transforms Fi{z) and F2{z), then Ri{a)R2{a) has a 
transform 




(6) 
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This will be true even if the transforms contain renormalon singularities (see 



Appendix). As noted by 't Hooft [11], the presence of such singularities in 



the factors leads to similar singularities in the product as a result of (|^). 
This aspect is discussed in more detail below. Indeed much of this paper 
can be considered a generalization of this observation about the convolution 
integral. 

The other straightforward operation is differentiation, since the transform 
of a"d"i?/da" is 



dz \ dz 

It is possible to remove powers of a from a^d^R/daP by using 

= F(0) + r e~''''F\z)dz. (8) 
a Jo 

The F{0) piece here corresponds to a delta function term in the transform 
which we will seek to avoid. Thus will only be used in this paper on 
functions for which -F(O) =0. 

The remaining operations 3-5 are all intrinsically difficult ones. By this 
we mean that it is unlikely that closed-form results will ever be found for their 
effect on transforms, in the way that and (0) exist for multiplication and 
differentiation. Even with a simple example like F{z) = 1/(1 — 2;), it appears 
that they can produce a transform that is much more complicated in its fine 
detail than the input one. Our simple and natural question does not have a 
simple and natural full answer. However the issue of practicality remains and 
for this a complete answer to the problem is not necessary. For the effective 
use of current and foreseeable QCD calculations it will be sufficient to extract 
only the more important effects that the use of these operations will have on 
the transforms. Furthermore, by concentrating on the main effects, a simple 
and natural picture does emerge. 

The general question of the effect of the operations 3-5 on Borel summable 
functions has been addressed before, in a pair of papers by Auberson and 



Mennessier [jT2[. Their main conclusion is that the functions thus produced 
are also Borel summable. This result is clearly a necessary starting point 
for our investigation in that it establishes that the transforms we will be 



trying to describe actually exist. In general [12| provides a firm foundation 



of rigour that will be relied upon implicitly throughout what follows. As 
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a consequence of the emphasis on extracting information that can be used 
in practical QCD contexts, the approach and language used here is however 
very different than in . 



3 The A-Expansion 

As a concrete example, consider operation 5, inverting a function R{a) to 
a{R). 

The central idea of our approach is to introduce the following split in the 
initial function 

R{a) = a + XR{a). (9) 

The parameter A is a bookkeeping device which will be set equal to one at 
the end. Thus 

oo 

Ria) = r^aJ^^^ (10) 

n=l 

is simply all of R{a) except for its lowest order term. Of course (|^) really 
defines a function i?(a. A) such that R{a) = R{a, 1). And if we are starting 
from a QCD observable, R{a) will be renormalisation scheme invariant, but 
R{a, A 7^ 1) will not be. But these nuances are unproblematical and the 
A-dependence of R{a, A) will be suppressed in the notation. 

It is also possible to think of a{R) as now having a A-dependence, though 
one more complicated than that in (0). Expand a{R) as an expansion in A. 



, , .da 
a{x) = X + A— 

OA 



+ 



A=o 2 dX^ 



+ ... (11) 

A=0 



The coefficients are functions of x and can be found by repeatedly differen- 
tiating (I) rewritten as 

X = a{x) + XR{a{x)) (12) 
with respect to A. Thus the first differentiation gives 

da ^ dR da , . 

and so 

% = -Ri^)- (14) 

OA A=0 
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The first few terms of (|TTD are found to be 

dR 

a(x)=x - XR(x) + X'^R(x) — 

dx 

This will be called the A-expansion of a{R). 

Given that all the other perturbative expansions in this paper are diver- 
gent, it is important to emphasize that the A-expansion is usually convergent 
as a series in A. A detailed discussion of this issue is deferred ||13[, but heuris- 



tic arguments can be offered for (^) being convergent. First, note that so 
far nothing has been assumed about the convergence properties of (p!0|). R{a) 
can be as favourable and well-behaved a function as required and ( [T5|) will 
still be derived. As a relation between R{x) and a{x) as functions, the struc- 
ture of the A-expansion is independent of their behaviour as expansions in 
X. In this light, there is no a priori reason to expect ( [T5| ) to be especially 
problematic for the cases of interest. Secondly, when R{x) is divergent, the 
resulting divergence of a{x) as a series in x can be thought of as having been 
absorbed into the R{x) that appear in (|15|) to leave a convergent expression. 
Thirdly, consider the function 

X{y) = X{y,x)^ y (16) 
R{x - y) 

and its inverse y{X) = y{X,x) for a fixed x. R{a) is expected to have one 
cut in the complex a-plane, along the negative real axis. Thus R{x — y) is 
analytic about y = for all a; > and X{y) will converge for \y\ < x. Hence 



its inverse, X{y), has a non-zero radius of convergence |Tj], but since 



a{x) = X — y{X,x), (17) 

so does (p!5|). 

Clearly none of these arguments is sufficient to prove that the A-expansion 



is convergent for A = 1, as required. Pending the detailed discussion [jT3[, the 
issue is set aside. For present purposes, we merely note that in the results 
presented here the expansions actually have infinite radii of convergence. 

Turning to the use of (^), for a given R{a) and hence R{a), setting A = 1 
in ( pTSD provides a systematic means of calculating a{R) from R{a). It is a 
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method of inverting functions. Its importance as such is that it serves to 
reduce the hard operation of inversion to a sequence (albeit infinite) of easy 
ones, namely differentiation and multiplication. Given the Borel transform of 
R{a), (ID - (ID can be used to calculate the Borel transform of any particular 
term in the A-expansion and then these contributions summed to obtain the 
Borel transform of a{R). 



4 Calculating the Transform 

The simplest case is where the transform of R{a) (and thus R{a)) has a single 
singularity at z = 1/a such that 

The generalization to multiple poles, as required for realistic QCD examples, 
will be straightforward. Neglecting numerical factors, the general term in 



(15) is 



A" . . . (19) 

gi + ...g„ = r2-l. (20) 

What is the transform of this, given an R implied by (|TB|) ? Using (0)-(§) 
this is easily found in principle, but for any particular F{z) that accords 
with (rS) the convolution integrals quickly become impossible to evaluate 



exactly. It becomes necessary to follow only the more important features of 
the transforms through the calculation. 

Note that if F{z) in (p!8| ) has another singularity at z = 1/a with 
P' = P — 1, then according to (§) the additional contribution to the factorial 
divergence of the coefficient r„ is suppressed by 1/n and so can be neglected 
at large orders. All singularities at 2; = 1/a with smaller (3 are similarily sub- 
leading in r„. The singularity with largest P at z = 1/a will be called the 
dominant singularity, the others sub-dominant. It is the dominant singular- 
ities (possibly at different z = 1/a) that will be most important in practice. 
Furthermore, for the operations 1-5 it turns out that sub-dominant terms in 
the initial transform only give rise to sub-dominant terms in the result. From 
now on the sub-dominant singularities will be neglected. 
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The features of the transforms that are to be tracked through the calcu- 
lation will be their behaviour close to z = 0, the positions of the poles and 
the dominant behaviour there. For any transform of interest, these can be 
summarised thus 



for an input transform like (p!8|). 

Note that all of these are also only singular at z = 1/a. Why is this ? 
Firstly, differentiation of R and the use of cannot cause the transform to 
become non-analytic at any other point. Sub-dominant terms also remain 
sub-dominant. Secondly, multiplication and the use of the convolution inte- 
gral have much the same effect. Consider where Fi{z) and F2{z) in (|^) 
are of the form (21). How does -^3(2;) behave ? For z < 1/a the integrand 



is finite and so is F3{z). But as 2; ^ 1/a the integrand begins to diverge at 
both ends of its interval and hence so can -^3(2) as z ^ 1/a. For z > 1/a, a 
principal value is taken in the integral where necessary (see Appendix) and 
^3(2;) is finite. The only point where -^3(2;) is non-analytic is 2; = 1/a. This 
is essentially the observation made by 't Hooft |]TT[. Furthermore, if -^1(2;) 
and -^2(2;) are as in (pl|), so is ^3(2;), but with 

= nil + + 1 (22) 

n ^ f -f3i + m2 + l if + ma + 1 < -/^a + mi + 1 

\ -(32 + mi + l if -/32 + mi + 1< -/3i + ma + 1 ^ ^ 

One has to know the behaviours near z = because when one of the trans- 
forms in the integrand of (^ is diverging to give the divergence in -^3(2;), the 
other transform's argument is tending towards 2 = 0. All this is neglecting 
the sub-dominant terms. However, such terms in either -^1(2;) or -Fa (2;) do 
not produce dominant terms in ^3(2;) as a result of this convolution; they 
are safely neglected. Since differentiation and multiplication are the only 
operations involved in finding ([T9|), its transform and any intermediate ones 
involved in finding it will indeed thus behave as (|2l|) . 

In addition, one can use (0), ( ^2]) and (^) to find the Ai, Pi and rrii of 
these transforms. Most of the terms ([19D turn out not to contribute to the 
dominant singularity in the transform of a{R). At each order in A there is 
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one term of the form 

y^P^n-l<^_^ (24) 

and it is only these that contribute to this singularity. Working through the 
(tedious) details and summing the A-expansion one finds that, for R{a) given 
by (|18|), the dominant part of the transform for a{R) is 

_ \^g-AT-i/a 

where ri is the one-loop coefficient of R{a). Finally, A is set equal to one. 

The main effect that inversion has had on the transform is thus to change 
the overall constant (the Ai in (PT|)). The /3j doesn't change; this need not 
be the case for other operations. Most importantly, the position of the sin- 
gularity hasn't changed. In retrospect this is an obvious consequence of the 
non-obvious fact that inversion can be reduced to multiplication and differ- 
entiation. 

Additional singularities at other positions do not change this basic pic- 
ture. If (|18D is generalised to 

i^(^)-E n_^"^ ^~V«™, (26) 

51) becomes 

m-Y. n z^l/a^. (27) 

The additional complications ( PU] ) introduces are all sub-dominant in (P7|). 

The idea of a split (^ in one of the functions leading to a A-expansion also 
renders operations 3 and 4 tractable. Indeed in these cases the A-expansion 
appears much more familiar. 

Ri{x + XMx)) = j:^R,ixr^ (28) 

n=0 

x + XR2{x) X X ) ' 

However since these operations involve two functions, the details and the 
results are contingent on the specifics of two input transforms and a discussion 
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is again deferred |T3|. A full set of results covers the situations one is liable 
to be confronted with in practice. 

A general feature is however clear. Because all these operations can be 
reduced to multiplication and differentiation, if the initial transforms have the 
renormalon structure predicted by QCD [Q 0, so do the transforms produced 
by the operations. The universality of that structure is further confirmed: 
all QCD observables look the same in the Borel plane. 
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Appendix 

Although it is a standard theorem that the Cauchy product of two Borel 



summable series is also Borel summable, the result for products required here 
is that the Borel transform of the product is given by the convolution (j^), 
even when singularities are present. Since we know of no proof of this in the 
literature, one is outlined here. Consider two transform integrals 

/(a) = r F{z)e-''''dz, g{a) = P H G{z)e-'/''dz, (30) 
Jo Jo 

where F{z) is analytic, but G{z) has a singularity at 2; = 77. The P indicates 
a principal value. The convolution integral will thus be 

H(z) = P r G{w)F{z - w)dw. (31) 
Jo 

Now 

h{a) = g{a)f{a) (32) 
= \wv( r ' G{z)e-^'^dz+ r G{z)e-''^dz\f{a) (33) 

= h{a)+\im dzG{z)e~'''' r dwF{w)e~'/''. (34) 

<:^OJo Jri-z-e 
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where _ 

h{a) = P / H{z)e-''''dz. (35) 

JO 

It is the final term in ( ^l]) that is at issue. The details of whether or not it 
tends to zero depend on the specific G{z). It suffices to consider whether 

e r ' G{z)dz (36) 

JQ 

vanishes as e ^ 0. If it does, then the natural generalisation of the convolu- 
tion theorem, namely that h{a) = h{a), is true. 
If, as if in QCD, 

Giz) ~ (r^ - z)-^, z ~ r/, (37) 

then ( p6D vanishes only if /3 < 2. However the naive principal value definition 

in (pOD only holds for (3 < 1 anyway. For /3 > 1, one can define the transform 

via _ 

g{a) = a-^'P / Giz)e-'/''dz (38) 
Jo 

for some n such that G{z) has a singularity with j3' = (3 — n < 1. In this 
paper, the transforms are implicitly defined like this. For clarity, the a~" 
factors are suppressed and /3 > 1 ia allowed, but with care it can be arranged 
such that (3 < 1 m all transforms. The convolution theorem thus holds for 
all products considered here. 

The generalization to both F and G having multiple renormalon singu- 
larities is straightforward. 

Principal values have been used to define the Borel integrals because the 
generalization of the theorem, (^) and (0), is then particularly natural. 
However it is more common in QCD to define Borel integrals using contours 
that detour around singularities. These versions have the disadvantage here 
that the contours required for the equivalents of (|3lD and (|35D are then not 
obvious. 
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